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Abstract 

The renormalizability of the self-avoiding manifold (SAM) Edwards model is established. 
We use a new short distance multilocal operator product expansion, which extends methods of 
local field theories to a large class of models with non-local singular interactions. This validates 
the direct renormalization method introduced before, as well as scaling laws. A new general 
hyperscaling relation is derived. Manifolds at the 0-point and long range Coulomb interactions 
are briefly discussed. 
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The statistical mechanics of fluctuating surfaces has attracted much attention in the 
recent years, with apphcations in many areas of physics, from string theories in high en- 
ergy physics to interface and membrane problems in soft condensed matter physics and 
biophysics In particular, tethered surfaces, which model polymerized flexible mem- 
branes, have unusual and interesting elastic properties. While these properties are now 
well understood theoretically for "phantom membranes" , that is when self-avoidance (SA) 
interactions are ignored, the consequence of incorporating SA constraints to describe real 
membranes is still an open problem. In practice, the search for a consistent theoretical 
treatment of SA interactions raises the fundamental question of applying renormalization 
group (RG) methods to extended objects, which is the issue addressed here. 

The theoretical study of SA polymerized membranes is centered around a model of 
tethered self-avoiding manifolds (SAM) [|[]^ directly inspired by the Edwards model for 
polymers The surfaces are generalized to intrinsically D-dimensional manifolds, rep- 
resenting L>-dimensional connected networks, whose nodes, labeled by internal continuous 
coordinates x G IR''^, are embedded in external d-dimensional space with position vector 
r(x) G IR*^. The associated continuum Hamiltonian 7i is 

^AbT = ^ Jd'^x (V,r(x))' + ^Jd''xJ rf^x' 5^(r(x) - r{x')) , (1) 

with an elastic Gaussian term and a self-avoidance two-body 5-potential with excluded 
volume parameter 6 > 0, non-local in "manifold space" IR^. 

A finite upper critical dimension (u.c.d.) d* for the SA interaction exists only for 
manifolds with a continuous internal dimension < D < 2. Phantom manifolds [b = 0) 
are crumpled with a finite Hausdorff dimension dn — 2D/{2 — D), and d* = 2dH- In 0-0,0 
an e-expansion about d* was performed via a direct renormalization (DR) method adapted 
from polymer theory [0. But many issues remain unanswered: The consistency of the DR 
method is proven only for D = 1 hy the famous mapping of (|l|) onto a (zero component) 
($^)^(r) field theory in external d-dimensional space 0. When D 1, model (P can no 
longer be mapped onto a local field theory, and the validity of RG methods and of scaling 
laws has been justified only at leading order through explicit partial resummations f^. 
The questions of a proper treatment for boundaries and of the value of the configuration 
exponent 7 [H are also open. 

In this Letter, we introduce a flexible formalism that allows us to prove the validity of 
the RG approach to self-avoiding manifolds, as well as to a larger class of manifold models 



1 



with non-local interactions. It broadly extends a recent work by the authors |^] for a simpler 



model [|T^, with a local singular interaction, of a phantom manifold interacting with a 



single impurity |]TT| . The present formalism is based on a new operator product expansion 



involving multilocal singular operators, and allows for a systematic analysis of the short 
distance ultraviolet (UV) singularities of the model. At the critical dimension d*, we can 
classify all the relevant operators and show that the model (|l]) is renormalizable to all orders 
by renormalizations (i) of the coupling b and (ii) of the position field r. As a consequence, 
we establish the validity of scaling laws for infinite membranes, as well as the existence of 
finite size scaling laws for finite membranes. The latter result ensures the consistency of 
the DR approach. A surprising result, which distinguishes manifolds with non-integer D 
from open linear polymers, is the absence of boundary operator renormalization, leading 
to the general hyperscaling relation 

1 = 1- i^d, (2) 

valid for finite SAM with D < 2, D ^ 1. Another surprise when considering SAM at the 
0-point is the appearance of a new relevant interaction term, which can supersede the 
usual three-body term. 

Perturbation Theory. For infinite SA manifolds, physical observables are expressible in 
terms of the P-point correlation functions, whose perturbative expansions are formally 



p 1 (^h\N r 2N P N 

(n^e''^-^^-)) = zj^^ wk y n^^^^^(n^^^'''^"\n '^^(r>2a)-?(x2.-i))>„ (3) 

The r.h.s. average (. . .)o is performed with respect to the ideal Gaussian manifold (6 = 0). 
The partition function Z in the denominator has a similar perturbative expansion in b, 
but with no external points. The product of 6 functions in (^) can be written in terms of 
exponential operators as 

N , r 2N d'^ki N- ^ 2N .r- . 

n 5^(r(x2„) -?(a:2a-i)) = / H ^ H CJkJ n e^'^^''^^^^ (4) 

with "dipolar constraints" Ca{ki} = (27r)^5'^(k2a-i + k2a) for momenta k^ G IR*^ (later 

called "charges") assigned to the points Xi. The correlation function (^) is defined as 

p ^ 

translationally invariant in external space, i.e. with the "neutrality rule" E qz = 0, a 

i=\ 



condition which is necessary when deahng with infinite membranes to avoid infrared (IR) 
singularities. The Gaussian average in is easily performed, using the identity 

^neii^.r(..)) = exp (-1 E k.k.G,,) (5) 

where Gij = —\xi — Xjp~^/((2 — D)Sd) is the massless propagator (Coulomb potential 
in D dimensions), with So = r{D/2) • Integration over the momenta then gives for the 
A^'th term of (|^) the manifold integral 

/ n d'^x, A-i exp ( - J E qzq^ ^] (6) 

J i=l \ ^ Z,m=l A J 

where A{xi} is the determinant associated with the quadratic form (now on IR) Q{ki} = 

2N 

S kikjGij restricted to the vector space defined by the N neutrality constraints Ca{ki}, 
k2a + ^2a-i = 0, and A;^ is a similar determinant involving also the external points zi 
and Zm 0. 

Note that a proper analytic continuation in D of (P) is insured from by the use 
of distance geometry, where the Euclidean measure over the Xi is understood as the cor- 
responding measure over the mutual squared distances distribution 
analytic in D. 

Singular Configurations. The integrand in (|^) is singular when the determinant 
A{xi} < 0. The associated quadratic form Q{ki}, restricted by the neutrality constraints 
Ca{ki}, is the electrostatic energy of a gas of charges ki located at Xi, and constrained to 
form neutral pairs a of charges (dipoles). For such a globally neutral gas, the Coulomb 
energy is minimal when the charge density is zero everywhere, i.e. when the non zero 
charges ki aggregate into neutral "atoms". When < D < 2, the corresponding minimal 
energy is furthermore zero, which implies that the quadratic form Q is non-negative and 
thus A > 0. Singular {xi} configurations, with A = 0, still exist when Q is degenerate, 
which happens when some dipoles are assembled in such a way that, with appropriate 
non-zero charges, they still can build neutral atoms. This requires some of the points Xi 
to coincide and the corresponding dipoles to form at least one closed loop (Fig. 1). This 
ensures that the only sources of divergences are short distance singularities, and extends 
the Schoenberg theorem used in 0. 

Multilocal Operator Product Expansion. A singular configuration can thus be viewed 
as a connected "molecule" , characterized by a set ^A of "atoms" p with assigned positions 
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Fig. 1: A general diagram with two external points and three internal dipoles 
(a); "molecules" describing singular conhgurations with one (b), two (c,d) and three (e) 
"atoms". (b,c,d) give UV divergences, (e) does not. 

Xp, and by a set C of links a between these atoms, representing the dipolar constraints Ca- 
For each p, we denote by Vp the set of charges i, at Xi, which build the atom p and define 
Ui = Xi — Xp for i G Vp. The short distance singularity of A~^/^ is analyzed by performing 
a small yi expansion of the product of the bilocal operators x') = 5'^{y{x) — r(x')) for 
the links a G £, in the Gaussian manifold theory (Eq. (0)). As will be shown below, this 
expansion around M. can be written as a multilocal operator product expansion (MOPE) 

n (p{x2a, X2a-l) = ^{Xp] C%...^ {Vi} (7) 

where the sum runs over all multilocal operators $ of the form: 

^{xp}= [d'^v n {■.{{V,y^5''{r-r{xp))} Apixp):] (8) 

Here Ap{xp) = A*^''p'^p^(Vx, r(xp)) is a local operator at point Xp, which is a product 
of x-derivatives of the field r, of degree Sp in r{xp) and degree Vp > Sp in V^;. (Vp)'^^ 
denotes a product of Qp derivatives with respect to r, acting on d'^{r — r{xp)). The symbol 
": :" denotes the normal product subtraction prescription at Xp (which, in a Gaussian 
average, amounts to setting to zero any derivative of the propagator Gij at coinciding 
points Xi = Xj = Xp). For Card(A^) = > 1, (||) describes the most general |7Vl|-body 
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contact interaction between the points Xp, with possible inserted local operators Ap{xp) at 
each point Xp. For \ = 1, it reduces to a local operator Ap{xp). 

The coefficient associated with the operator $ in the MOPE, ^{j/i}, can be written 
as an integral over the momenta k^: 



^ ^ Jo.ec peM l^ieVp ieVp 

(9) 

where C^p{yi, k^} is a monomial associated with the operator Ap, of similar global degree 
Tp in the y^, and Sp in the kj. The product H' is over all constraints a E C but one. 

The MOPE (0) follows from the expression (^) in terms of free field exponentials plus 
constraints. For each p, we use the general small yi local operator product identity 



i,3evp (10) 



When expanded in the y^, the normal product :( )\xi=xp'- in (|T0|) gives a sum E C'^jyi, k^} 
:A(xp)e^'^p''p : (denoting kp = E k^ and Yp = r(xp)) which generates the local operators 

A{xp) and the monomials C"^ of (P) and (|[). We insert the identity 1 = / (i^kp5'^(kp — 

S ki) in (^) for each atom p E M., rewrite one of the dipolar constraints as a global 

ieVp 

neutrality constraint 5^( E kp) on the kp, and expand each 5'^(kp — E k^) in powers 

peM ieVp 

of kp. Finally by integrating over the kp, the constraint d'^{ E kp) builds the multilocal 

peM 

|7V/i|-body operator ^{xp} and we obtain the MOPE (|^), (D and (|). 

Power Counting and Renormalization. The MOPE (|^) allows us to determine those 
singular configurations which give rise to actual UV divergences in the manifold integral 
(^. Indeed, given a singular configuration Ai and integrating over the domain where the 
relative positions yi = Xi — Xp are of order we can use the MOPE of (§) to obtain an 

expansion of the integrand in (^) in powers of p. Each coefficient ^ gives a contribution 
of order p'^, with degree to given by power counting 



L'{2|/:|-|A^|} + cfz^o{|A^|-|/:|-l}+ E {MQp-Sp)+rp} (11) 

pEM 



with t'o = (2 — -D)/2 < 1 and Vp > Sp. Whenever w < 0, a UV divergence occurs, as a factor 
multiplying the insertion of the corresponding operator $. At the upper critical dimension 
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d* = 2D/uo, u becomes independent of the number \C\ of dipoles, and is equal to the 
canonical dimension a;$ of f lid^x^ in the Gaussian theory. Three relevant operators, 

M 

with ujci, < and such that the corresponding coefficient does not vanish by symmetry, 
are found by simple inspection. Two of these operators are marginal (c<;$ = 0): (i) the 
two-body SA interaction term d'^{rp — Tp') itself, obtained through singular configurations 
with = 2 atoms (and with q = r = s = for p and p'), (ii) the one-body local elastic 
term :(Vrp)^:, obtained for = 1 {q = 0, r = s= 2). The third operator is relevant with 
u!<i, = —D, and is just the identity operator 1 obtained when \A4\ = 1 {q = r = s = 0). It 
gives "free energy" divergences proportional to the manifold volume, which cancel out in 
IR finite observables (|^). 

The above analysis deals with superficial UV divergences only. A complete analysis 
of the general UV singularities associated with successive contractions toward "nested" 
singular configurations can be performed, using the techniques of and the fact that an 
iteration of the MOPE only generates multilocal operators of the type (^. The results are: 
(i) that the observables are UV finite for d < d*{D), and are meromorphic functions 
in d with poles at d = d*, (ii) that a renormalization operation, similar to the subtraction 
operation of can be achieved to remove these poles, (iii) that this operation amounts 

to a renormalization of the Hamiltonian (|1]) . More explicitly, the renormalized correlation 

p 

functions ( II e^^''"^'-^''')„ have a finite perturbative expansion in the renormalized coupling 

7=1 

6r, when ^- ■ ■'j^ is the average w.r.t the renormalized Hamiltonian 
nn/ksT = f / rf''x(V.r^(x))' + ^6^^^Z, J rf^x J rf^x' 5^(r^(x) - r^(x')) • (12) 

is a renormalization (internal) momentum scale, e = 2D — duo, Zi,{bii) and Z^b^) are 
respectively the coupling constant and the field renormalization factors, singular at e = 0. 
At first order, we find by explicitly calculating cjp^'^'^ and C^^ that Z = 1 + 6^-^ ^^^^-* , 
Zb = 1 + r(2D/(2-gj) ^ B = i(47r)-^^^+^(2 - D)-^+^ . For the quantities 

which are not IR finite, which we discuss later, an additive counterterm proportional to 
the volume of the manifold (corresponding to the relevant identity operator 1) is also 
necessary. 

Expressing the observables of the SAM model (P in terms of renormalized variables 

one can derive in the standard way RG equations involving 
Wilson's functions VF(6r) = ^^-^bJ^,\^, z^(6r) = 1^0 — ^^j.-^\nZ\^. A non-trivial IR fixed 
point &R oc e is found for e > 0. It governs the large distance behavior of the SA infinite 
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manifold, which obeys scahng laws characterized by the exponent defined for instance 
through the 2-point function ((r(x) — r(0))^) oc \x\'^" . The value obtained in this approach, 
f = z^(6^), corroborates that obtained in |@-||,|| at first order in e. 

Finite Size Scaling and Direct Renormalization. The DR formalism requires one 
to consider finite manifolds with "internal volume" V , and to express scaling functions in 
terms of a dimensionless second virial coefficient g = —K~^Z2^cl{,Z\)^i where 2i{V) and 
2^2, c(^) are respectively the one- and two membrane (connected) partition functions, and 
Rq is the radius of gyration. 

When dealing with a finite closed manifold (for instance the iD-dimensional sphere 
So 0), characterized by its (in general curved) internal metric, the massless propagator 
Gij gets modified. Nevertheless, from (^0]) and the short distance expansion of Gij in a 
general metric [|T2| , one can show that the short distance MOPE (^ remains valid, provided 
that the sum is extended to include multilocal operators $ still of the form (^, but with 
local operators A{x) involving also the Riemann curvature tensor and its derivatives, with 
appropriate coefficients ^. A crucial point is that in the MOPE the dependence on 
the geometry of the manifold (size, curvature,. . .) is encoded only in the expectation values 
^. . . $ . . of the multilocal operators $, while the short distance behavior [yi 0) of 
coefficients ^{yi} is independent of the geometry. Thus, at (i*, UV divergences still 
come with insertions of relevant multilocal operators with uj^ < 0. 

When < D < 2, none of the new operators involving the curvature is found to 
be relevant by power counting. Therefore, the infinite membrane counterterms Z and 
still renormalize the finite membrane theory. Since, as for finite size scaling [jl^, the 
manifold size is not renormalized, arguments parallel to those of for polymers can be 
used to justify the DR formalism. Indeed, the second virial coefficient g{h, V) (as well 
as any dimensionless scaling function) must be UV finite when expressed as a function 
Snibin ^/U^) of the renormalized coupling (and of /u). As a consequence, (i) the scaling 
functions are finite when expressed in terms of g and obey RG equations. The existence of 
a non-trivial IR fixed point 6* for e > implies that (ii) in the large volume limit V oo, 
g tends toward a finite limit g* = g^ib'^) (independent of Vn^), and so do all scaling 
functions. Points (i) and (ii) are the essence of DR. 

Hyperscaling. As mentioned above, the renormalization of partition functions for a finite 
SAM requires an additional counterterm (shift of the free energy) proportional to the 
manifold volume V. A consequence of the absence of other geometry dependent relevant 
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operators when < D < 2 is the general hyperscahng law valid for closed SAM, and 
relating the configuration exponent 7, defined by 

Z^{V) = J V[r\ 5'^{r{0)) q-^/^bT _ ^ ^^3^ 

to the exponent v. Indeed, from (|1^), once the free energy divergent term has been 
subtracted, Zi is simply multiplicatively renormalized as Ziih, V) = Z~'^l'^ Z^{h^, V/j^). 
This vahdates the scaling hypothesis that Zi ~ I?!""' ~ y-i^d/D^ leads directly to (|2|). 

For open SAM with free boundaries, and when 1 < D < 2, the boundary operator 
r d^~^x 1 becomes relevant. Since it is simply a geometrical quantity, it cannot modify 

boundary 

the renormalizations of r and b. Furthermore, it is marginally relevant only for D = 1 
and therefore, as long as D ^ 1 the scaling laws and the hyperscahng relation (^) remain 
valid. Only at -D = 1, the corresponding end-point counterterm enters the multiplicative 
renormalization of Zi, and 7 becomes an independent exponent, with an extra contribution 
from the two end-points. 

Eq. (|^) has been checked explicitly at order e for the sphere So and the torus Tp. Pre- 
vious calculations which yield 7 = 1 for non-integer D, did not involve the physical 
massless propagator [jl2| (valid for a finite manifold with Neumann boundary conditions), 
used here. 

When D > 2, if the small e RG picture remains valid, i.e. if the large distance 
properties of SAM are governed by the IR fixed point 6* , operators involving curvature 
become relevant, and (|^) is not expected to be valid, even for closed manifolds. 

0-point and long-range interactions. The above formalism is directly applicable to a 
large class of manifold models where the interaction can be expressed in term of free field 
exponentials with suitable neutrality constraints Ca{ki}. Examples of such interactions are 
the n-body contact potentials but also the two-body long-range Coulomb potential l/|r — 
r'!*^"^, which can be represented by modified dipolar constraints C{ki} = |k|~^5'^(k -f k'). 
For all these models, the MOPE involves the same multilocal operators as in (^, with 
modified coefficients (still given by (H), but with new constraints Ca)- 

As an application, we may ask for the most relevant short-range interaction describing 
a polymerized membrane at the 0-point, i.e. when the two-body term b in (P vanishes. 
It is either the usual three-body contact potential, with u.c.d. ^3 = 3D/{2 — -D), as 
for ordinary polymers, or the two-body singular potential Ap5'^(r — r') with u.c.d. ^2 — 
2(3-D — 2)/(2 — D), which indeed is the most relevant one when D > 4/3. 
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Finally, the absence of long-range potentials in the MOPE shows that long-range 
interactions are not renormalized. For instance, when considering charged polymerized 
membranes with a two-body Coulomb potential, the only (marginally) relevant operator 
at the u.c.d. is the local operator :(Vr)^:, indicating that only r is renormalized. As a 
consequence, it is easy to show that v = 2D/ {d — 2) exactly in this case. 
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